Langevin equation approach to granular flows in narrow pipes 
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The flow of granular material through a rough narrow pipe is described by the Langevin equation 
formalism. The stochastic force is caused by irregular interaction between the wall and the granular 
particles. In correspondence with experimental observations we find clogging and density waves in 
the flowing material. 
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When granular material flows through a narrow ver- 
tical pipe one observes recurrent clogging and density 
waves . This effect is well known to physicists and 
engineers, usually it is undesirable and causes technologi- 
cal problems, e.g. in chemical engineering. Density waves 
play a major role in the behavior of granular materials 
and have been investigated by many authors using vari- 
ous methods: 

Ristow and Herrmann Q] reproduced the density fluc- 
tuations in an out-flowing hopper by molecular dynamics 
that have been observed experimentally before (e.g. ||) 
using the X-ray technique. Baxter and Behringer [g] sim- 
ulated the flow using an cellular automaton. Peng and 
Herrmann J5J studied a Lattice Gas Automaton Q for 
the flow of granular material. Using phenomenologically 
plausible rules for the interaction of particles and of par- 
ticles with the wall they could reproduce density fluc- 
tuations which spectrum obeys a power law. Lee and 
Leibig ||] applied the kinetic wave approach JlIJ to the 
flow of granular particles through a pipe. They treated 
initial random density fluctuations as a set of distinct 
homogeneous density regions and considered the motion 
of the interfaces between them. Finally they could show, 
that the evolution of such a simple model leads to the 
formation of clusters with a high density contrast. 

The aim of the present paper is to provide a model 
for the flow of granular material in a vertical narrow 
pipe using the Langevin-equation approach of stochastic 
forces. Such a approach was also used by Mehta et al. p"l| ] 
to describe the relaxational behavior of a granular pile 
submitted to vibration and was proved to be suitable. 
We will show that our model is able to reproduce the 
experimental observations QJ^]. The simulation of the 
resulting density equation requires much less computa- 
tional effort than the direct simulation using molecular 
dynamics IpqJlSfl. We discuss the instability of the ho- 
mogeneous flow in the hydrodynamic approximation and 
provide critical values for the occurrence of clogging and 
density waves. 

When sand flows through a narrow pipe we assume 
that there is a permanent random interaction of the sand 
particles with the wall of the pipe. The equations of mo- 



tion for a single particle which is subjected to gravity g 
in positive x-direction and which does not interact with 
other particles in the low density regime read 



(la) 
(lb) 



The friction 7 and the Langevin fluctuation term results 
from the interaction of the grain with the wall. (Although 
our model does not include the interaction of the grains 
with the air inside the pipe, in a very simple approxima- 
tion one can assume that the fluctuation term accounts 
for this interaction too.) For the stochastic force we as- 
sume Gaussian white noise ((£,i(t)£,j(t + T)) = 5ijS(T)). 
The path of a particle is scattered independently at dif- 
ferent places during its motion downwards which is de- 
scribed by independent impacts in time. Hence, after the 
relaxation time m/7 the velocities of the particles obey 
a Maxwellian distribution with mean value v° = mg/'y. 

We apply the collision integral proposed by Prigogine 
and Herman |b| for the description of particle interac- 
tions. In their investigation it was intended to model ve- 
hicular traffic flow and it has been pointed out by several 
authors (e.g. [[liPl) that traffic flow on one-lane high- 
ways reveals striking similarities to granular flow in a 
pipe. When a fast moving particle i collides with a slower 
one j after the collision both grains move with the lower 
velocity Vj. 
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During the impact the momentum balance of the par- 
ticles is not conserved. Due to the strong interaction 
between the particle and the wall we suppose that the 
lost part of the momentum will be taken over by the wall 
which is assumed to be coupled with a reservoir. Intro- 
ducing the effective cross section C which depends on 
the geometry of the pipe and the particles the collision 
integral reads 



dP 
~dt 



P(x,v,t) P(x,v',t) (v'-v) dv' 



C P(x,v,t) n(x,t) (u(x,t) - v) , (2) 
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where P(x,v,t) is the one-particle probability density in 
phase space. The particle density n(x, t) and the mean 
velocity u(x, t) at position x and time t are given by 



n(x, t) 
u(x, t) = 



P(x, v, t) dv 



1 



n(x, t) J _ 
Hence we find the kinetic equation 
dP 



v P(x, v, t) dv 



(3a) 
(3b) 



dt dx ^ dv 



P 



ej d 2 P 



(g - —v) 

+ CP(x, v, t)n(x, t)(u(x, t) — v) 



ra 2 dv 2 

and in a homogeneous stationary flow (n = n , u = u° 

m („. „.o\ 2 



(4) 



n° exp 



2 k B T 



2ir k B T a 
The mean velocity 

rag Ck B T°n° 
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(5) 



(6) 



depends on the density and the mean square displace- 
ment of the velocity: 



T(x,t) = 



(v — u(x,t)) P(x,v,t)dv. (7) 



k B n(x, t) 
In the homogeneous case we find 

T = T° = — 
k B 

Therefore the homogeneous flux through the pipe 
/ rag C k B T° 



f = n°u° 
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(8) 



(9) 



shows two distinct regimes: a low density regime due to 
a high particle velocity where only very few collisions oc- 
cur and a high density regime due to low particle velocity 
caused by dissipative impacts of particles. 

When we assume local equilibrium, i.e. n° — » n(x,t), 



i(x,t), and T° — * T(x,t), we rewrite eq. 



dP d 
dt dx 



yP] 



d_ 

dv 



F(n,T) 7 



v P 



e 7 d 2 P 



m 2 dv 2 
(10) 

Similar as in the Vlassov formulation flij] the force 

F(n,T) = mg~C k B T(x,t)n(x,t) (11) 

is determined to be self-consistent in its dependence on 
density and granular temperature of the material. With 
the effective force F(n, T) acting on the particles at a 



given location x and a given time t we find the corre- 
sponding Langevin equation for the motion of the parti- 
cles which are subjected to gravity and impacts of other 
grains 

Xi = Vi (12) 
mvi = —yvi + F (n(x l ,t) , T (x t , t)) + y/2ej &(t) . 



Eqs. (|11|,|12|) determine an effective simulation algorithm 
(for details see [pill ). 
From eq. (|l0|) we derive the hydrodynamic equations 
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(13a) 
(13b) 
(13c) 



The first two terms on the right hand side of the heat 
balance equation (13c) describe the heat exchange be- 
tween the granular material and the wall, whereas the 
last term leads to an effective volume viscosity. 

In the approximation of quick temperature and ve- 
locity relaxation one finds for the high damping limit 
(7 — > 00) the Burgers equation 



dt 7 dx 7 dx 2 



(14) 



In this limit there are no self sustained inhomogeneous 
solutions of eq. ( |l4| ) (see |l6[]). For finite damping, 
7 < 00, however, as shown below the homogeneous solu- 
tion n = n°, U = u°, T = T° of the eqs. ( p"3| ) becomes 
unstable when the average density approaches a critical 
value n cr . In the context of clustering instabilities in 
dissipative gases Goldhirsch and Zanetti M] argued sim- 
ilarly: when the pressure in a dense region decreases due 
to dissipation, the resulting pressure gradient leads to 
further increase of the density which finally results in a 
granular cluster. 

In the following we want to discuss the stability anal- 
ysis of the hydrodynamic equations ( 1^) |Q . Obviously, 
there is a homogeneous solution (pip ) for a given homo- 
geneous density n° . We disturb the homogeneous state 
in eqs. (|l3| ) with wave-like perturbations (Sn ~ du ~ 
ST ~ exp(— at + ikx)), drop the quadratic terms and 
get a eigenvalue-problem for a(k). If the real part of 
a(k) is negative, Re(a(fc)) < 0, fluctuations can grow 
and the initially homogeneous state becomes unstable. 
The transition occurs at the critical particle density n cr 
where Re(a(fc)) = 0. Because of the assumed peri- 
odic boundary conditions the wavenumber k is discrete: 
k = ^jji (i — ±1, ±2, ■ • •). Fig. |l| shows the critical den- 
sity n cr over the dimensionless mode number i. Obvi- 
ously in particular short length perturbations are able to 
destabilize the homogeneous flow in the granular system. 



2 



This stands in strong contrast to results found for hydro- 
dynamic formulations for traffic p3| , ^9|j20| and granular 
flows |2^] , where the long range fluctuations are the crit- 
ical ones. Our results are not surprising if one imagines 
that a local large gradient of the velocities will lead to a 
high collision rate at this place. For sufficient high den- 
sity this process leads to clusters with high local density 
and small average velocity. 

For large wave numbers we get a low limiting critical 
density given by 



lim 



'(ft) 



7 



7 



7 



1 



(15) 



where Lb = V e m /7 is the braking distance, i.e. the 
length a particle is damped out after an impact. It deter- 
mines the length scale which characterizes our granular 
system and its critical behavior. In contrast, the critical 
behavior of the traffic flow models proposed in p3| , ^9| , p0[ 
depends on the length L of the entire (periodic) system 
too since the critical fluctuations are long range ones. 

To check the analytic results the time-discretized 
Langevin equations 

Xi(t + M)=Xi(t)+Vi(t)At (16a) 



V2e7 At 



GRND (16b) 



have been solved numerically. GRND is a Gaussian ran- 
dom number with standard deviation equal unity. Using 
the parameters m = 7.4 10 -7 kg, 7 = 7 10 -6 kg/sec, 
e = 2.0 KT 8 Nm, C = 6.4 1(T 3 , g = 9.81 m/sec 2 , and 
At = 0.01 sec we calculated the time dependent density, 
velocity and granular temperature from the trajectories 
of the Brownian particles. The given parameters have 
been determined experimentally |l8|. We assume peri- 
odic boundary conditions and homogeneous initial con- 
ditions. Inserting the given parameters into eq. ( )l5| ) we 
find n cr = 12000/m. 

Fig. shows the velocity distributions 



w(v,t) = / P(x,v,t)dx 
Jo 



(17) 



of a stable (undercritical) (n° = 11000/m < n cr ) and 
an unstable system (n° = 14000/to > n cr ). In the un- 
dercritical case we find a stable homogeneous flow with 
Gaussian velocity distribution, while in the latter case in- 
homogeneities due to random fluctuations increase with 
time and the velocity distribution w(v, t) is no longer 
Gaussian. In our opinion there are at least two distinct 
velocity distributions in the system: at regions of low 
density we find high average grain velocity and at high 
density regions the grains move with low average veloc- 
ity. The overlayed curves in the distribution plots (fig. |^) 



show support of the assumption of a bimodal velocity 
distribution. 

Fig. H shows snapshots of the particle density of the un- 
stable system where gravity acts in positive a;-direction. 
We eventually observe the formation of two moving clus- 
ters originating from random inhomogeneities in the early 
state of the system-evolution. In dependence from the 
initial conditions we found also configurations with one 
or three moving clusters. In correspondance with the ex- 
periment and MD simulations Q we observed coexisting 
clusters moving either in positive or negative direction. 
Note that the slope at the left hand side of the density 
wave in fig. |^ is very steep. Here the collision rate is very 
high due to the large velocity gradient between the parti- 
cles which are involved in the clusters and the free falling 
ones (fig. ||). The particle velocity at the right hand 
side of the clusters is much lower. There the grains' ve- 
locity slowly increases under the influence of gravity and 
hence the high density area, i.e. the cluster, dissolves (see 
also pH|). These processes lead to the obvious hump-like 
shape of the clusters in fig. ||. Contrary to the hump- 
like solutions of the Burgers-equation [|16| the widths of 
the clusters remains invariant when they move through 
the pipe. Fig. [| shows the described sharp decrease and 
slowly increase of the particle velocitie s. A s one can de- 
rive from the hydrodynamic equations (13c ) the high neg- 
ative velocity gradient at the front (left hand side) of a 
cluster leads to an increase of the granular temperature, 
whereas the small positive velocity gradient inside and 
at the backside of the cluster results in a smaller granu- 
lar temperature compared with outside the clog. Fig. |] 
(lower part) illustrates this plausible behavior. 

We investigated the granular flow trough a vertical nar- 
row pipe using a simple model consisting of Brownian 
particles with collision interaction. We could show that 
there is a critical value for the particle density which de- 
cides whether a initially homogeneous flow remains sta- 
ble. The model is valid in the limit of pairwise particle in- 
teraction. This precondition is assumed to be fulfilled for 
the case of moderate particle density and low pipe width. 
Our model also defines an efficient simulation method for 
granular flows. Applying this algorithm to low- and high- 
density pipe flow we found that in both cases the numer- 
ical results agree with the theoretical prediction. The 
numerical results for the spacial particle density, the av- 
erage velocity, and the granular temperature agree with 
the hydrodynamic description. 

We thank S. Esipov, S. Savage, J. Schuchhardt, and 
H. J. Herrman for useful discussion. 
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FIG. 1. 
number i. 



The critical density n cr over the mode 
The parameters are m = 7.4 1CT 7 kg, 



7 = 7 l(r b kg/sec, e = 2.0 1(T S Nm, C = 6.4 V0T s , 
g = 9.81 m/sec 2 . For diflerent parameters the curve 
changes, however, its qualitative shape remains con- 
served. 
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FIG. 2. Velocity distribution found by simulations of 
the stable homogeneous and the unstable system (solid 
lines) and the corresponding analytical results (dashed 
lines) . 
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FIG. 3. The particle density of the unstable system at 
times t = Asec and t = 500sec. Since the initial (homoge- 
neous) density is overcritical the inhomogeneities increase 
with time. 
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FIG. 4. The mean velocities (top) and the mean 
square displacement of the velocity (bottom) of the un- 
stable system at time t = 500sec. 
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